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Summary  of  Project  Goals. 

Based  on  the  known  principles  that  govern  protein  folding  of  polypeptide 
chains  into  antiparallel  p-sheet  structures,  a  smaH  protein  molecule,  p-bellin,  has 
been  designed  by  Drs.  Jane  and  David  Richardson.  Amino  add  sequertoes  of 
several  p-^in  variants  have  been  synthesized  by  Dr.  Bruce  Erickson.  In  the 
current  project,  computer  methods  were  used  to  evaluate  the  structure  of  p-beltin 
models  and,  more  generally,  to  answer  fundamental  questions  about  the  relations 
among  molecular  geomei^.  atomic  potential  energy,  and  static  stability  of  p- 
sheeted  proteins.  Using  CONGEN,  a  program  for  conformatiored  search,  macro- 
molecular  energy,  minimization  and  dynamics  calculations,  the  stereochemistry, 
molecular  surface,  p-sheet  twist  and  properties  of  p-bellin  were  areilyzed  and 
compared  with  those  of  other  typical  p-sheeted  proteins  such  as  immunoglobu¬ 
lins,  prealbumin  or  superoxide  dismutase.  The  project  was  particularly  directed 
towards  improving  the  unsatisfactory  features  of  solvent-accessible  surface  of  the 
P-bellin  molecule.  Quenched  dynamics  simulations  were  used  to  improve  on  sur¬ 
face  smoothness  aruj  side  chain  interactions. 


Summary  of  Accomplishments 

(1)  Preliminary  characterization  of  p-beilin  model. 

Coordinates  of  the  model  of  p-bellln  were  obtained  from  Jane  Richardson 
and  Richard  Feldmann.  Potential  energy  of  the  model  was  examined  by  the  pro¬ 
gram  CHARMM,  version  16.  The  energy  of  the  model  (P-bellin  dimer  with  566 
atoms)  was  comparable  to  that  of  an  unrefined  crystall^raphic  structure  (=  1 
kJ/atom).  The  geometry  of  the  model  was  not  optimal  but  readily  improved  upon 
energy  minimzation. 

A  disulfide-bonded  version  of  p-beltin  was  consbucted  next  by  substituting 
side  chains  of  Val  21  for  cysteines  (sulphur  atoms  were  taken  to  be  identical  to 
the  CGI  carbons  of  the  Val  residue).  This  construction  did  not  introduce  any  large 
strain  into  the  structure.  After  200  cycles  of  tiie  ABNR  minimization  protocol,  the 
total  energy  of  the  structure  improved  to  =  -1700  kJ.  Although  the  Cp-S-S-CP  tor¬ 
sional  angle  was  not  optimal,  it  seemed  clear  that  manual  adjustment  of  tNs  tor¬ 
sion  to,  e.g.,  90"  followed  by  additional  energy  minimization  would  produce  a 
flawless  disulfide. 

The  Chothia-Janin-Teiler  formula  for  the  exposed  surface  of  a  single-domain 
protein  gives  the  expected  accessible  surface  for  p-bellin  39.96  nm^.  The  sur¬ 
face  obtained  using  the  Lee  &  Richards  (1971)  algorithm  on  the  (unminimized) 
model  was  45.13  nm^,  i.e.  13%  larger.  This  difference  may  be  indicative  of  an 
incorrect  fold.  On  the  other  harxl,  the  ratio  of  nonpolar/poiar  surface  side  chain 
atoms  for  P-beiiin  was  1 .2  while  it  was  about  1 .7  for  native  proteins  and  above  2.0 
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for  inoorrsctiy  fokjed  models. 

The  last  property  of  the  ^-bellin  model  ttiat  was  systematically  examined  was 
the  p-sheet  twist  Haight  lines  were  fitted  into  the  ^strand  tMCkbones  by  a 
least-square  procedure  and  the  angles  between  successive  p-strand  axes  were 
measured.  All  the  strand-strartd  angles  were  exceptionaity  low  (4^;  2.4*;  5.6*)  and 
remained  so  after  minimization.  The  unusual  flatness  of  p-beliin  barrel  can  be 
visualized  by  least-squares  fitting  of  a  twisted  hyperboloid  into  the  ^-strands. 
Such  surfaces,  when  fitted  to  common  ^-barrels,  have  semiaxes  of  eilipticai 
cross-sections  of  about  0.9  nm  and  0.6  nm.  resftoctiveiy,  with  helical  pit<^  in 
between  8-9  nm.  The  ^-bellin  hyperboloid,  on  the  other  hand,  was  characterized 
by  semiaxes  lengths  2.5  nm  and  0.5  nm,  and  its  helicai  pitch  was  very  high  (61 .3 
nm)  indicating  virtually  no  twist. 

(2)  Molecular  representation  of  the  BAB  linker. 

One  of  the  more  successful  p-bellin  preparations  contained  a  bivalent  linker, 
BAB,  which  tethers  the  two  C-termini  of  the  ^-bellin  dimer.  The  linker  is  a  com¬ 
plex  organic  molecule  combining  aryl  derivatives  with  motifs  of  peptide  chemistry. 
In  order  to  represent  correctly  the  structure  and  stability  of  the  BAB-IInked  p-beliin 
dimer,  a  "Residue  Topology  RIe*  and  "Parameter  RIe"  of  the  BAB  molecule  were 
prepared,  the  BAB  molecule  representation  was  generated  in  the  CONGEN  pro¬ 
gram,  and  its  stereochemistry  and  energetics  was  analyzed.  Dr.  R.  E.  Bruccoleri 
partidpated  in  BAB  molecule  construction. 

(3)  Analysis  of  backbone  shapes,  backbone  twist  and  side  chain  interactions 
in  p-sh^ed  protein  domains. 

One  of  the  approaches  to  the  problem  of  protein  folding  has  been  to  study 
protein  shapes.  Single  domains  built  of  antiparallel  P-etrands  were  altematively 
described  as  twisted  prisms  (C.  Chothia),  etliptical  barrels  (J.  Richardson)  and 
twisted  hyperboloids  (J.  Novotny,  R.  Bruccoleri  &  J.  Newell).  We  have  completed 
a  mathematical  least-squares-fitting  study  of  higher-power  curved  geometrical 
surfaces  into  polypeptide  badtoones  of  8  different  ^-barrels.  Current  models  of 
^-barrels  involve  the  supposition  of  a  regular  average  right-harxJed  twist  of  the  ^ 
strands;  we  have  found  that  large  deviations  from  the  average  p-sband  twist 
values  occur  in  positions  of  polar  side  chains.  In  superoxide  dismutase,  the 
exceptional  backbone  twist  values  (exceeding  6  »  420*  ±20*)  are  so  common  that 
the  supposition  of  a  regular  backbone  twist  has  to  be  abandoned.  A  profile  of  side 
chain  energies  Is  dominated  by  large  fluctuations,  showing  the  striking  anisotropy 
of  atomic  force  distribution  throughout  the  structure.  The  regular  right-handed 
twist  seems  to  be  characteristic  of  the  norvpolar,  spherically  symmetric  van  der 
Waals  forces  among  side  chain  atoms,  whereas  hydrogen  bonding,  with  its  strict 
directionality,  tends  to  discrupt  this  pattern.  The  picture  of  ^-barrel  energetics 
which  emerged  from  the  analysis  is  an  unexpected  one.  Regular  p-barreis  can  be 
constructed  from  coiled  arto  irregular  backbones,  and  deviations  from  the  ideal, 
regularly  twisted  p-sheet  geometry  are  quite  common.  The  importance  of 
geometrical  models,  such  as  the  twisted  h^>erboloid,  thus  becomes  apparent:  an 
abstract  mathematical  model  facilitates  an^sis  of  P-barrei  properties  arto  hek>s 
to  bridge  the  gap  between  the  geometrical  and  physical  attributes  of  the  struc¬ 
tures. 
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(4)  Electrically  charged  molecular  surface  of  p-beUln. 

Following  a  suggestion  of  Dr.  Jane  Richardson,  electrically  charged  solvent 
accessK)le  surfaces  were  computed  in  a  series  of  proteins,  arul  compared  with 
the  percentage  of  electrically  charged  surfaces  in  various  ^-bellins.  The  results 
showed  that  the  soluble  preens  to  have  at  least  10%  of  their  solvent-exposed 
surface  charged.  The  original  ^-bellin  model  had  7.3%  charged  surface,  whereas 
the  BAB-IInked  p-bellin,  with  dramatically  improved  solubility,  has  1 1 .2%  charged 
surface. 

(5)  Quenched  molecular  dynamics  simulations. 

Quenched  dynamics  consists  of  molecular  dynamics  simulation  at  a 
predetermined  temperature,  followed  by  a  gradual  cooling  of  the  structure.  It  is 
presently  the  most  powerful  energy  minimization  protocol  available  for  computer 
manipulation  of  pro^  structures. 

In  the  first  run.  the  heating  phase  of  the  simulation  involved  3  ps  of  10*  K 
temperature  increments  every  100  fs  spanning  ttie  temperature  range  from  0*  K 
to  300*  K.  During  each  temp^ature  increase  the  velodties  of  atoms  were  reas¬ 
signed  according  to  a  Gaussian  distrtoution.  Ten  ps  of  simulation  were  then  car¬ 
ried  out,  followed  by  3  ps  of  cooling;  the  cooling  represented  a  reversal  of  the 
heating  protocol,  /.e.  10*  K  temperature  decrements  every  100  fs  from  300*  K  to 
0*  K.  The  second  run  involved  5  ps  heating  phase,  from  0*  K  to  500*  K.  10  ps  of 
dynamical  simulation  at  500*  K,  and  5  ps  of  cooling  down  to  0*  K.  Figure  1  shows 
the  time  course  of  total  energy  (light  line)  and  temperature  (heavy  line)  changee 
in  the  simulations.  Several  interasttog  points  emerged  from  the  analy^  of  the 
completed  dynamics  runs. 

(I).  Relatively  large  movements  were  apperent  throughout  the  whole  structure, 
including  the  BAB  linker.  At  higher  temperature  (500*  K)  the  N-terminal  ^ 
strand  showed  tendency  to  untold  away  from  the  rest  of  the  structure. 

(H).  Significant  rearrangements  were  visible  at  the  tight  turns  which  Include  pro¬ 
line  reeiduee.  Hers,  the  backbone  twist  angle  6  changed  quite  dramatically, 
oompared  to  the  starting  ^-belll^  structure. 

(W).  Perhaps  most  bnportantly,  the  relatively  fiat  p-beNin  sheet  acquired  a  distinct 
right-harrded  sheet  twist  on  quenched  dynamics.  This  was  in  contrast  to  the 
M-handed  sheet  twist  found  in  ali^the  naturally  occuring  p-sheets.  It  seems 
that  the  driving  force  of  sheet  twis^  originated  from  the  proHne-contalning 
turns.  Thus,  analysis  of  energetics  and  geometry  of  the  quenched  ^-belilns 
concentrated  on  the  possible  role  of  ^tums  In  P-sheet  architecture. 

(6)  Molecular  dynamics  simulations  to  mplore  p-tum  conformations. 

In  order  to  get  insight  into  the  role  of  side  chains  in  p-tum  conformation, 
dynamicai  simulations  were  run  not  only  on  the  original  p-bellin  model  but  also  on 
its  sequence  variants  where  proline  residues  were  replaced  with  either  glycines  of 
D-prolines.  The  results  of  the  simulations  demonstrated  convincingly  the  relation¬ 
ship  between  p-tum  twist,  side  chain  chirality,  and  the  hydrogen  bonding  pattern 
inside  the  turn.  The  turns  7-8-9-10  and  15-16-17-18  of  the  original  p-bellin  coordi¬ 
nates  had  the  canonical  hydrogen  bonding  pattern  (see  also  Fig.  2)  whereas  the 
turn  23-24-25-26  had  an  abnormal  hydrogen  bond  between  the  residues  24  and 
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26,  more  reminfsoent  of  Y*tum8  (see  also  Fig.  2).  TNs  difference  seemed  to 
predetermine  all  the  24-25  turns  to  become  type  I’  (lefthanded,  see  Fig.  3)  while 
the  other  turns  acquire  either  the  righthanded  type  I  (Fig.  4)  or  the  type  I’  twist 

The  greater  backbone  flexibility  of  8.16,24-Qly-p-beilin  was  reflected  in  larger 
fluctuation  of  turn  twist  values.  Nevertheless,  the  trends  observed  in  p-beilin 
were  reproduced  in  the  Gly-p-beilin.  That  is,  the  type  I  twist  was  stabilized  by  the 
side  chain  carboxylate  hydrogen  bonding,  and  t^  I*  turns  tended  to  dispense 
with  the  unusual  Phe-Pro  N  -»  O  hydrogen  bond  as  their  twist  became  more 
negative. 

Introduction  of  the  D-proline  into  p^jellin  turns  was  accompanied  by  negative 
(type  I'-iike)  twists.  Turns  which,  upon  dynamics,  maintained  the  backbone- 
backbone  hydrogen  bonds,  or  acquired  the  backbone-side  chain  carboxylate 
hydrogen  bonds,  tended  to  acquire  more  positive  twists.  The  type  I*  twist  seemed 
to  be  stable  even  when  backbone-backbone  hycfrogen  bonds  vanished. 

(7)  Molecular  dynamics  simulations  with  impiidt  solvent  effects. 

During  the  last  part  of  the  funding  period,  the  quenched  molecular  dynamics 
simulations  were  repeated  under  conditions  that  simulated  more  realisticaily  sol¬ 
vent  eftects.  That  is,  (i)  the  non-bonded  interactions  were  evaluated  to  infinity;  (il) 
the  AMBER/CHARMM  19  parametrization,  with  larger  van  der  Waals  radii  and 
larger  partial  atomic  charges  (simulating  p^arization  effects)  was  employed;  (lU) 
and  the  electrostatic  potential  function  had  the  form  £J«(0i02)/t16icR^), 
where  the  "dielectric  constanT,  0  •  4R,  better  approximates  the  electrostatic 
screening  effects  of  bulk  solvent  This  set  of  simuiations  produced  more  con¬ 
sistent  results,  without  abnormal  p-tum  twists  seen  previously,  arto  without  side 
chairvmain  chain  hydrogen  bonding  in  toe  turns.  As  a  rule,  p-Mlin,  contakiing  L- 
Pro,  tended  to  unfold.  For  example,  the  root-mean-square  deviation  between  the 
starting,  energy-minimized  structure  and  the  final  "quenched”  structure  was  0.5 
nm  after  10  ps  at  500"  K,  compared  to  the  r.m.s.  deviation  of  only  0.3  nm  for  the 
Qly-  and  D-Pro-p-bellins  under  the  same  conditions.  The  significcmtiy  higher  sta¬ 
bility  of  the  Qly-  arto  D-Pro  p-bellins  correlated  with  the  opposite  chiral  terKferKies 
of  their  respective  P-tums.  Thus,  L-Pro  imposed  right-harKied  twists  on  turns 
(423"  average  value),  while  Gly  or  D-Pro  imposed  left-handed  twists  (averages  of 
-1"  and  -6",  respectiveiy).  This  set  of  molecular  dynamics  simuiations  appeared 
to  establish  a  direct  role  of  p-tum  twist  (and,  impii^y,  side-chain  chirality  of  the 
turns)  in  p-sheet  stability. 

(8)  Criteria  that  discriminate  between  native  proteins  and  incorrectly  folded 
mod^. 

Various  theoretical  concepts  such  as  free  energy  potentials,  electrostatic 
interaction  potentials,  atomic  packing,  solvent-exposed  surface,  surface  charge 
distribution  etc.,  were  tested  for  their  ability  to  discriminate  between  native  pro¬ 
teins  and  misfolded  protein  models.  Misfolded  models  were  constructed  by  intro¬ 
ducing  incorrect  side  chains  onto  polypeptide  backbones:  side  chains  of  the  a- 
heiical  hemerythrin  were  modeled  on  toe  Mheeted  backbone  of  immunoglobulin 
VL  domain  while  those  of  the  VL  domain  ware  similarly  modeled  on  the  hem- 
erytorin  backbone.  CONGEN,  a  conformational  space  sampling  program,  was 
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us«d  to  construct  the  side  chains.  CapabHity  of  the  conformationai  search  pro¬ 
cedure  to  reproduce  native  conformations  was  gauged  first  by  rebuilding  (the 
correct)  side  chains  in  hemerythrin  and  the  VL^  domain:  constructs  with  r.m.s.  j 

differences  from  the  X-ray  side  chains  2.2  -  2.4  A  were  produced,  and  many  cal¬ 
culated  conformations  matched  the  native  ones  quite  well.  Incorrectfy  folded 
models  were  then  constructed  using  the  same  conformational  protocol  applied  to 
incorrect  amino  add  sequences.  All  CONQEN  constructs,  both  corrs^  and 
incorrectly  folded,  were  characterized  by  exceptionaily  small  molecular  surfaces  * 

and  low  potential  energies.  Surface  chai^  density,  atomic  packing  and  Coulomb 
formula-based  electrostatic  interactions  of  the  misfolded  structures  and  the 
correctly  folded  protekw  were  similar,  and  therefore  of  little  Interest  for  diagnosing 
incorrect  folds.  The  following  criteria  dearly  favored  the  native  structures  over  the 
misfolded  ones:  (1)  solvent-exposed  side  chain  non-pdar  surface,  (10  number  of  I 

buried  ionizable  groups,  arfo  (ilO  empirical  free  energy  functions  that  incorporate 
solvent  effects. 
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Giy-^>b«l]ln  minimized  end  after  SOO  K  dynamics,  r.m.s.  3.1  A 
Implicit  solrent  by  AMBER/CHARlIlf  19  parametrlzation  &  D^Ar*  elec. 
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